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INTRODUCTION 

A.  History 

As  knowledge  of  the  ear's  response  to  intense  sounds  has 
developed  over  the  years,  the  desire  for  a  good  theoretical 
understanding  of  the  processes  involved  has  been  amplified 
by  a  practical  concern  for  predicting  or  ameliorating  the  ef¬ 
fect  of  intense  noise  exposures.  A  substantial  portion  of  the 
world’s  population  is  exposed  to  intense  impulses  in  a  mili¬ 
tary  setting;  but  sport  shooting  and  other  less  bellicose  civil¬ 
ian  activities  also  result  in  intense  noise  exposures  to  an  even 
larger  number.  And  there  is  no  question  that  even  brief  expo¬ 
sure  to  intense  impulses  can  result  in  permanent  loss  of  hear¬ 
ing.  Therefore,  the  need  for  valid  information  with  respect  to 
impulse  noise  hazard  is  critical  because  the  population  at 
risk  is  large  and  the  possibility  of  damage  is  real. 

There  are  a  number  of  criteria  for  impulse  noise  expo¬ 
sure  in  use  in  the  world  {China,  France,  Germany,  Great 
Britain,  The  Netherlands,  and  the  United  States  of  America 
all  have  their  own  criteria  [CHABA,  1968;  Cheng  et  al. , 
1987;  MIL  STD-1474  (B),  1979;  Ministry  of  Defense  (Brit¬ 
ish),  1982;  Ministry  of  Defense  (French),  1982;  Pfander, 
1975;  Smoorenburg,  1980]}.  However,  all  the  criteria  have 
been  derived  essentially  empirically,  with  little  theoretical 
basis.  So  long  as  we  have  only  a  rudimentary  understanding 
of  the  mechanisms  acting  in  the  ear  at  high  sound  pressures, 
attempts  to  find  the  right  combination  of  measures  to  use  in 
rating  hazard  will  be  severely  hampered.  Furthermore,  the 
latest  data  on  hazard  suggest  that  the  applicability  of  all  the 
criteria  is  seriously  in  question  in  that  they  over  rate  the 
hazard  from  low-frequency  impulses  (Dancer  et  al.,  1983; 
Price  et  al.,  1989). 

B.  Applied  issues 

In  order  for  a  criterion  for  noise  exposure  to  be  useful,  it 
should  rate  hazard  accurately.  In  addition,  it  is  an  advantage 
if  it  is  formulated  on  sound  theoretical  bases  in  acoustics, 
physiology,  and  psychophysics.  The  issues  of  accuracy  and 


form  are  separable.  An  empirically  derived  criterion’s  pre¬ 
dictions  could  correlate  highly  with  hazard  for  the  specific 
type  of  sound  used  in  its  derivation  (and  therefore  be  accu¬ 
rate  for  that  type  of  impulse);  but  it  may  not  be  generalizable 
to  other  types  of  stimulation  and  may  provide  no  physical 
insight  into  the  processes  that  are  responsible  for  damage. 
On  the  other  hand,  if  a  criterion  was  theoretically  based, 
then  in  addition  to  rating  hazard  for  specific  sounds,  it  would 
possess  the  important  benefit  of  generalizability  to  the  effects 
of  sounds  not  used  in  its  derivation,  Furthermore,  it  could 
have  heuristic  value  in  suggesting  designs  or  procedures  that 
would  reduce  or  eliminate  hazard.  Clearly,  a  theoretically 
based  criterion  would  possess  real  advantages. 

An  adequate  understanding  of  the  ear’s  response  to  im¬ 
pulse  noise  should  provide  answers  to  three  perplexing  ques¬ 
tions  that  challenge  our  current  concepts.  First,  why  are 
low-frequency  impulses  (Friedlander  waveforms  with  long 
A  durations)  not  as  dangerous  as  impulses  in  the  midrange? 
For  equal  peak  pressures,  the  low-frequency  impulses  con¬ 
tain  much  more  total  energy,  and  they  both  contain  the  same 
energy  in  the  midrange;  yet  as  we  have  seen  earlier,  the  low- 
frequency  impulses  are  less  damaging.  Second,  why  is  it  that 
audiometrically  measured  or  histologically  established 
losses  occur  in  the  midrange,  even  for  impulses  that  have 
spectral  peaks  near  100  Hz?  Third,  why  are  losses  not  much 
greater  than  they  actually  are  in  practice?  That  may  seem 
like  a  perverse  question;  but  consider  the  great  energy  pres¬ 
ent  in  some  impulse  exposures.  Just  one  impulse  from  a  large 
cannon  could  contain  as  much  A-weighted  energy  as  a  full 
year  of  exposure  to  a  hazardous  industrial  noise.  Often  firers 
wear  no  hearing  protection  while  many  rounds  are  fired. 
Why  then  are  the  losses  not  much  greater  than  they  actually 
are? 

C.  Operating  assumptions 

Given  the  requirement  for  an  accurate  understanding  of 
the  basis  of  the  problem  as  well  as  the  apparent  inadequacy 
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of  existing  methods  for  rating  hazard,  some  years  ago  we 
initiated  a  research  program  intended  to  address  the  situa¬ 
tion.  Presently,  our  working  hypotheses  are  that  within  the 
organ  of  Corti  there  is  a  critical  displacement  or  stress  at 
which  the  mechanism  of  loss  becomes  essentially  mechani¬ 
cal  in  nature.  Once  that  level  is  reached,  additional  stimula¬ 
tion  results  in  a  rapid  growth  of  damage.  Furthermore,  the 
free-field  sound  pressure  producing  the  critical  stress  is 
largely  a  function  of  the  spectral  tuning  of  the  external  and 
middle  ears,  which  effectively  form  a  bandpass  filter  (Price, 
1986b). 

We  hasten  to  point  out  that  the  models  to  be  presented 
in  this  paper  are  not  definitive;  but  are  still  in  active  develop¬ 
ment.  However,  they  do  provide  answers  to  the  foregoing 
questions,  and  also  provide  a  ready  basis  for  the  develop¬ 
ment  of  research  hypotheses. 

I.  THE  MODELS 
A.  General  considerations 

For  the  purposes  of  this  paper,  we  will  try  to  present 
enough  technical  detail  with  respect  to  the  modeling  to  make 
the  implications  of  a  model  intelligible.  An  earlier  prelimi¬ 
nary  report  on  the  model  has  been  published  (Kalb  and 
Price,  1987)  and  full  account  of  the  modeling  is  presently  in 
preparation. 

Several  general  considerations  guided  the  modeling. 
First,  we  wanted  to  be  able  to  relate  the  various  elements  of 
a  model  to  specific  physiological  entities  and  to  give  them 
physically  realistic  values,  to  allow  for  maximum  physical 
insight.  Second,  to  insure  maximum  utility,  we  wanted  the 
model  to  be  calculable  with  modest  computer  resources  (a 
PC-level  computer).  Third,  the  main  focus  of  the  model  was 
on  the  behavior  of  the  ear  at  high  sound  intensities  where  loss 
was  thought  to  be  a  function  of  mechanical  stress  (higher 
than  120  or  130  dB  SPL).  Last,  the  model  had  the  goal  of 
beginning  with  a  pressure  history  in  the  free  field  and  passing 
the  energy  through  the  outer,  middle,  and  inner  ears  to  pro¬ 
vide  an  estimate  of  displacement  within  the  inner  ear.  In 
order  to  meet  these  last  goals,  we  had  to  integrate  models  of 
the  head,  outer  ear,  middle  ear,  and  inner  ear  and  allow  for 
nonlinear  elements  in  them.  We  gratefully  acknowledge  the 
efforts  of  many  investigators  who  have  produced  models  of 
individual  sections  or  components  in  this  chain.  As  a  conse¬ 
quence  of  their  efforts,  we  have  primarily  been  concerned 
with  the  selection  and  integration  of  components  into  a  uni¬ 
fied  system. 


B.  The  models 

The  details  of  the  models  presented  are  appropriate  for 
the  cat  ear.  This  choice  was  made  because  the  physiological 
and  acoustic  values  were  best  known  for  the  cat  and  also 
because  much  of  the  hearing  loss  data  to  be  explained  were 
produced  with  the  cat  ear.  The  ears  of  other  mammals  are 
highly  similar  in  structure  and  the  principles  used  in  the 
modeling  should  transfer  with  only  modest  adaptations  to  fit 
their  anatomic  details. 

The  ear  was  modeled  with  electroacoustic  elements. 
Mathematically,  development  proceeded  in  two  alternate 
forms,  one  in  the  time  domain  and  one  in  the  frequency  do¬ 
main  (Kalb  and  Price,  1987).  The  time  domain  model  can 
be  characterized  as  equations  of  motion  formulated  as  a  sys¬ 
tem  of  coupled  ordinary  differential  equations  that  allow 
both  linear  and  nonlinear  calculations  in  the  time  domain. 
The  frequency  domain  version  can  of  course  only  be  applied 
to  linear  processes.  Where  the  nonlinearity  of  the  stapes  was 
treated  (discussed  later),  we  performed  a  nonlinear  trans¬ 
formation  of  the  data  and  then  proceeded  with  the  remaining 
elements  of  the  model.  Both  models  share  the  same  set  of 
variables  and  values  so  that  their  performance  is  compara¬ 
ble.  The  two  models  are  complimentary  in  that  they  each  are 
structured  to  allow  efficient  computation  of  particular  pa¬ 
rameters.  The  circuit  diagram  of  the  external  and  middle  ear 
portions  of  the  models  appears  in  Fig.  1. 

1.  From  the  free  field  to  the  cochlea 

The  head  was  modeled  as  a  spherical  baffle  around 
which  the  sound  wave  diffracts  and  in  which  the  ear  canal  is 
located.  The  acoustic  field  up  to  the  ear  canal  entrance  was 
approximated  by  a  simple  network  proposed  by  Bauer 
( 1967)  and  the  Wiener  el  at.  ( 1965)  two  tube  model  of  the 
ear  canal  follows  it.  The  middle  ear  elements  are  based  on  a 
model  by  Zwislocki  (1962)  with  additional  impedance  ele¬ 
ments  given  by  Lynch  etal.  ( 1982).  These  elements  account 
for  the  fact  that  the  ear  drum  does  not  move  as  a  unit,  for 
losses  in  the  ligaments  joining  the  ossicles,  for  the  areal  ratios 
of  the  ear  drum  and  stapes  and  the  lever  ratios  of  the  ossicles, 
and  for  the  resonances  of  the  volumes  included  in  the  middle 
ear. 

We  can  gain  confidence  in  the  performance  of  the  inte¬ 
grated  model  by  comparing  values  calculated  from  the  mod¬ 
el  with  acoustical  quantities  that  have  in  fact  been  measured. 
If  for  the  moment  we  allow  the  cochlea  and  its  loads  to  be 
represented  as  linear  mass  and  resistance  elements,  then  the 


FIG.  1.  Circuit  diagram  of  model. 
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model’s  calculations  can  be  compared  with  the  available 
data.  The  various  transfer  functions  for  the  ear  appear  in 
Figs.  2-5.  Figure  2  is  the  transfer  function  from  free-field 
pressure  to  ear  drum  pressure  (data  points  from  Wiener  et 
ai,  1965);  Fig.  3  is  the  transfer  function  of  the  middle  ear 
(data  from  Guinan  and  Peake,  1967);  Fig.  4  shows  the  im¬ 
pedance  of  the  ear  (data  from  Tonndorf  and  Khanna,  1967, 
and  Moller,  1963);  and  Fig.  5  shows  intracochlear  pressure 
as  a  function  of  pressure  at  the  drum  (data  from  Nedzel- 
nitzky,  1980).  The  agreement  between  the  models  and  the 
data  is  good  except  for  frequencies  above  10.0  kHz  in  Fig.  2 
and  above  5.0  kHz  in  Fig.  4.  We  believe  the  lack  of  agree¬ 
ment  in  these  areas  is  due  primarily  to  measurement  difficul¬ 
ties  inherent  at  high  frequencies  (noted  by  the  original  au¬ 
thors).  In  general,  we  believe  that  the  integrated  model  does 
a  creditable  job  of  reproducing  the  quantities  that  can  in  fact 
be  measured  at  low  intensities. 

2.  The  nonlinear  stapes 

The  stapes  is  commonly  modeled  as  a  linear  element. 
However,  at  very  high  intensities  (above  130  dB ) ,  the  mid¬ 
dle  ear  becomes  observably  nonlinear  (Guinan  and  Peake, 
1967;  von  Bekesy,  1936).  We  hypothesize  that  the  primary 
reason  for  the  nonlinearity  is  the  annular  ligament,  which 
holds  the  stapes  in  place  at  the  entrance  to  the  cochlea,  and 
because  of  its  mechanical  structure  can  exert  a  controlling 
influence  on  cochlear  input  (Price,  1974).  A  preliminary 
report  has  been  given  (Price  and  Kalb,  1986),andadetailed 
description  of  the  mathematical  model  of  the  stapes  suspen¬ 
sion  and  a  full  consideration  of  its  effects  are  being  written 
for  publication. 

Strictly  speaking,  the  assertion  that  the  annular  liga¬ 
ment  is  responsible  for  the  nonlinearity  goes  somewhat  be¬ 
yond  the  available  data.  However,  the  ligament  is  a  likely 
candidate  for  the  distinction.  In  comparison  with  the  rest  of 
the  middle  ear,  the  ligament  is  capable  of  withstanding  im¬ 
mense  mechanical  forces.  Of  course,  it  is  possible  for  other 
parts  of  the  middle  ear  to  become  nonlinear.  For  example,  in 


FIG.  2.  Comparison  of  calculated  values  (solid  lines)  with  measurement  of 
the  transfer  function  of  the  outer  ear  by  Wiener  el  al.  ( 1965 ). 


FIG.  3.  Comparison  of  calculated  values  (solid  lines)  with  measurement  of 
the  transfer  function  of  the  middle  ear  by  Guinan  and  Peake  ( 1967). 


some  animals,  the  malleo-incudal  joint  is  capable  of  slippage 
(even  if  it  is  almost  never  observed),  the  incudostapedial 
joint  could  pull  apart  or  slide,  and  the  ear  drum  itself  is  clear¬ 
ly  capable  of  greater  outward  than  inward  movement.  The 
list  could  go  on;  but  we  believe  that  these  nonlinearities  are 
likely  to  be  preceded  by  the  nonlinearity  at  the  stapes,  and  as 
a  consequence,  their  effect,  so  far  as  the  cochlea  is  concerned, 
is  overridden  by  the  limitation  imposed  by  the  annular  liga¬ 
ment.  In  any  event,  whatever  the  source  of  the  nonlinearity, 
the  result  is  reduced  efficiency  in  energy  transmission  into 
the  cochlea. 

The  annular  ligament  model  conception  in  mechanical 
form  is  shown  in  Fig.  6.  The  stapes  is  viewed  as  suspended  at 
the  edges  of  its  footplate  by  hinged  springs,  a  configuration 
that  naturally  allows  for  high  compliance  when  displace¬ 
ments  are  small  (as  is  true  for  the  normal  range  of  sounds  to 
which  ears  are  exposed);  but  which,  because  of  the  geome¬ 
try,  quickly  becomes  stiffer  as  displacement  increases.  The 
values  for  the  spring  suspension  were  chosen  to  fit  data  from 


FIG.  4.  Comparison  of  calculated  values  (solid  lines)  with  measurement  of 
the  impedance  of  the  ear.  Data  points  from  Tonndorf  and  Khanna  (1967) 
and  Moller  ( 1963). 
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FREQUENCV  IN  KHZ 

HG.  5.  Comparison  of  calculated  sallies  (solid  lines)  of  intracochlear 
pressure  as  a  function  of  pressure  at  l he  drum  with  measurements  of  Ned- 
/clnil/ky  ( |9X0). 


Lynch  el  al.  (1982)  on  static  and  dynamic  stiffnesses  for  the 
stapes.  In  addition,  as  displacements  rise  we  make  the  stiff¬ 
ness  and  damping  increase  proportionally.  The  increases  in 
damping  are  just  sufficient  to  keep  the  middle  ear  from  'ring¬ 
ing'  when  it  is  driven  hard.  The  full  effect  of  the  suspension 
geometry  and  spring  characteristics  can  be  seen  in  Fig.  7, 
which  portrays  the  restoring  force  for  any  given  displace¬ 
ment.  As  indicated  in  the  figure,  the  stapes  displacement 
) caches  an  asymptote  at  20 gm  ( "breaking  displacement"). 

In  Fig.  8,  the  model’s  output  is  compared  with  actual 
measurements  of  middle  ear  displacement  made  by  Guinan 
and  Peake  ( 1967)  in  the  anesthetized  cat.  It  can  be  seen  in 
this  figure  that  the  model  fits  the  available  data  and  limits 
maximum  displacement  to  about  20 gm  in  either  direction. 

3.  The  cochlea 

The  cochlea  is  modeled  as  a  two-chambered,  fluid-filled 
box  with  rigid  side  walls.  The  partition  between  the  cham¬ 
bers  is  rigid,  except  for  a  tapered  basilar  membrane  ( BM) 


UNLOADED  L0A0ED  OEflECTEO  BREAKING 

ANNULAR  LIGAMENT  MODEL  STATES 
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FIG.  7  Diagram  showing  the  calculated  restoring  force  on  the  stapes  as  a 
function  of  displacement 


that  becomes  progressively  more  compliant  toward  the  apex 
( mechanical  stiffness,  mass,  and  viscous  losses  changing  ex¬ 
ponentially ).  In  addition,  the  cross-sectional  area  of  the 
cochlea  decreases  exponentially  from  base  to  apex. 

For  rapid  computation,  we  used  a  WKB  approximation 
to  calculate  transfer  functions  to  each  of  50  points  along  the 
BM.  By  using  forward  and  inverse  FFTs,  we  calculated  a 
history  of  BM  motion  along  its  length.  In  the  differential 
equation  model,  the  cochlea  is  discretized  into  512-coupled 
resonant  elements  whose  characteristics  change  with  loca¬ 
tion.  The  values  chosen  for  the  mechanical  properties  pro¬ 
duce  the  envelopes  of  traveling  wave  displacements  that 
show  a  changing  Q  and  an  increase  in  amplitude  with  in¬ 
creasing  frequency  as  seen  by  von  Bekesy  ( llM9).  We  believe 
that  these  envelopes  are  characteristic  of  the  BM  when  driv¬ 
en  at  high  amplitudes  where  the  sharp  tuning  produced  by 
active  processes  that  operate  at  low  intensities  cannot  have 
an  appreciable  effect. 


:0a  :  is  120  130  H0  150  160  170  isa 

RMS  EARDRUM  SOUND  PRESSURE  LEVEL  ( DB  RE  0.0000  D1NES-CM-2) 


FIG.  6.  Diagram  of  the  suspension  of  the  stapes  in  the  annular  ligament  FIG  8,  Comparison  of  calculated  stapes  displacement  with  measurements 
model  showing  the  various  states  it  can  assume.  by  Guinan  and  Peake  ( 1967). 
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4.  Hazard  calculation 

The  mathematical  models  were  first  used  to  calculate 
the  envelopes  of  BM  displacement  (Kalband  Price,  1987)  as 
a  means  of  assessing  relative  hazard  of  different  impulses. 
But  the  BM  undergoes  many  oscillations  in  response  to  any 
but  the  simplest  stimulus;  therefore  to  account  for  the  effects 
of  the  repeated  flexing  of  the  BM  during  an  exposure,  we 
added  an  additional  calculation.  As  a  first  approximation, 
we  assume  that  loss  processes  operating  within  the  cochlea 
are  akin  to  those  that  are  generally  operative  when  materials 
are  repeatedly  stressed  and  finally  fatigue  and  fail  (Price, 
1983a).  As  Broch  (1979)  suggested,  for  materials  in  gen¬ 
eral,  the  number  of  cycles  to  failure  is  a  function  of  the  strain 
raised  to  a  power,  usually  greater  than  two,  i.e.,  N  xSh  —  K, 
where  A"  is  a  constant  and  b  is  the  power  to  which  the  strain 
is  raised.  It  is  apparent  that  this  formula  could  be  integrated 
with  the  model’s  calculation  of  BM  displacement  to  yield  an 
estimate  of  hazard.  Therefore,  we  calculated  a  value  (strain) 
for  the  peak  of  the  upward  displacement  of  each  of  the  50 
segments  of  the  BM  model  as  it  was  exposed  to  a  given  stimu¬ 
lus,  the  value  was  raised  to  some  appropriate  power  (2.5  was 
used  in  the  calculations  reported  here)  and  the  effect 
summed. 

An  additional  factor  needs  to  be  noted  regarding  the 
calculation  of  strain.  The  BM’s  width  changes  about  3  fold 
from  base  to  apex;  therefore  it  follows  that  the  amount  of 
strain  represented  by  any  given  displacement  must  also  vary 
as  a  function  of  the  linear  dimension  of  the  member  being 
displaced.  To  account  for  this,  in  calculating  hazard  we 
weighted  displacements  of  the  BM  progressively  less  with 
more  apical  locations  in  accordance  with  the  change  in  BM 
width.  The  end  result  of  this  exercise  was  a  number  for  each 
of  the  50  segments  of  the  BM  that  integrated  both  the  peak 
amplitude  of  its  displacement  with  the  number  of  occur¬ 
rences.  As  this  is  written,  we  have  only  begun  to  explore  the 
hazard  calculation.  The  units  of  hazard  are  arbitrary,  al¬ 
though  we  presume  that  the  calculated  value  represents  at 
least  an  ordinal  scale  of  hazard. 

In  calculating  hazard,  only  the  upward  deflections  of 
the  BM  were  included.  This  was  done  because  such  displace¬ 
ments  had  earlier  appeared  to  correlate  best  with  the  damage 
( Kalb  and  Price,  1987)  and  have  the  additional  advantage  of 
making  physiological  sense  as  possible  damage  producers. 
Briefly,  upward  movements  of  the  basilar  membrane  pro¬ 
duce  predominantly  stretching  forces  at  the  level  of  the  tec¬ 
torial  membrane  and  reticular  lamina.  Given  that  tissue 
tends  to  fail  in  tension,  upward  displacements  could  well  be 
especially  hazardous.  Conversely,  compressive  movements 
(downward  displacements),  would  be  expected  to  be  less 
damaging.  A  recent  conception  of  the  micromechanics  of 
transduction  in  hair  cells  ( Pickles  et  al.,  1984)  fits  this  line  of 
argument.  Upward  displacement  of  the  basilar  membrane  is 
associated  with  forces  that  cause  the  stereocilia  to  deflect 
toward  the  longer  stereocilia,  stretching  the  tip  links,  de¬ 
creasing  the  membrane  resistance  and  exciting  the  hair  cell. 
Movements  in  the  opposite  direction  simply  allow  the  links 
to  go  slack.  If  this  is  in  fact  the  case,  then  stretching  move¬ 
ments  (upward  displacements)  could  easily  result  in  dam¬ 
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age  and  compressive  movements  (downward  displace¬ 
ments)  could  be  relatively  less  hazardous. 

II.  APPLICATION  OF  THE  MODELS 
A.  Effect  of  the  middle  and  external  ear 

The  influence  of  the  conductive  properties  of  the  middle 
ear  on  an  impulse  is  immense.  The  bandpass  filtering  effect  is 
well  known  and  the  proposed  peak  limiting  add  to  it.  The 
effect  of  peak  limiting  can  be  appreciated  by  considering  the 
displacement  that  would  have  taken  place,  had  the  middle 
ear  been  linear.  In  response  to  impulses  in  the  crew  areas  of 
large  caliber  weapons,  calculations  indicate  that  the  stapes 
would  try  to  move  2000  or  more  microns.  With  peak  limiting 
present  (as  is  the  case  in  the  real  ear),  the  displacements  are 
limited  to  less  than  20  fim,  a  reduction  of  more  than  two 
orders  of  magnitude! 

The  effect  of  the  middle  and  external  ears  on  the  energy 
that  enters  the  cochlea  is  portrayed  in  Fig.  9.  In  this  some¬ 
what  complicated  figure  the  energy  in  1/3  octave  bands  en¬ 
tering  the  cochlea  has  been  calculated  for  Friedlander  wave¬ 
forms  with  A-durations  characteristic  of  a  rifle  and  a 
howitzer,  each  at  three  intensities.  These  waveforms  are 
close  simulations  of  weapons  impulses;  but  do  not  have  the 
“hash”  on  them  associated  with  measures  of  real  impulses. 
Two  major  points  can  be  made  from  this  calculation.  First, 
even  though  the  spectral  peaks  of  these  impulses  in  air  would 
have  been  at  140  and  800  Hz,  respectively,  at  the  input  to  the 
cochlea  they  all  have  their  greatest  energy  at  about  3000  Hz. 
The  ear  acts  like  a  fairly  sharply  tuned  bandpass  filter;  there¬ 
fore  all  impulses  naturally  tend  to  have  their  greatest  effects 
in  the  ear’s  mid-range.  Second,  note  that  with  increases  in 
intensity,  the  energy  increases  nonlinearly  for  both  impulses. 
However,  because  of  the  peak  limiting,  the  howitzer  type 
impulse  at  180-dB  peak  pressure  actually  has  less  energy  in 
the  midrange  than  the  160-dB  rifle-type  impulse.  This  pro¬ 
cess  could  explain  much  of  why  the  low-frequency  impulses 
are  less  hazardous  than  previously  thought. 
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FIG.  9,  Calculated  1/3  octave  energy  entering  the  cochlea  for  Friedlander 
waveforms  simulating  rifle  and  howitzer  impulses  at  three  intensities. 
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B.  Effects  for  the  intact  ear 

Given  that  we  are  interested  in  modeling  the  ear’s  re¬ 
sponse  to  intense  impulses  and  the  attendant  hazard,  it  is 
instructive  to  explore  what  happens  when  the  model  is  driv¬ 
en  with  acoustic  waveforms  that  have  been  used  in  experi¬ 
ments  with  real  ears.  These  calculations  can  be  seen  in  Figs. 
10-15.  In  figures  where  the  hazard  function  is  displayed,  the 
maximum  on  the  vertical  scale  for  the  hazard  function  was 
arbitrarily  set  at  10,000  “units”  in  order  to  allow  at  least 
relative  comparisons  between  impulses. 

1.  The  166-dB  howitzer  impulse 

In  Fig.  10,  we  see  the  response  of  the  ear  to  a  166-dB 
howitzer  impulse  used  in  a  previous  study  (Price,  1986a). 
The  top  panel  is  the  pressure  history  of  the  impulse,  the  mid¬ 
dle  panel  is  the  calculated  stapes  displacement,  and  the  bot¬ 
tom  panel  is  the  calculated  BM  displacement  envelope  and 
the  hazard  function.  It  is  apparent  that  while  the  calculated 
stapes  displacement  preserves  the  general  form  of  the  im¬ 
pulse,  the  effect  of  the  clipping  is  visible.  All  the  extremes  of 
pressure  appear  compressed  in  the  stapes  displacement  and 
the  oscillations  near  ambient  pressure  are  proportionately 
larger,  as  would  be  expected. 

The  envelope  of  BM  displacements  is  relatively  broad 
and  shows  one  peak  in  midcochlea  and  a  second  one  near  the 
apex.  When  the  evolution  of  the  BM  displacements  resulting 
from  the  stapes  motion  was  examined,  then  a  surprising  ef- 


FIG  iO.  166-dB  howitzer  impulse  (Price,  1986a),  The  upper  panel  is  the 
pressure  history,  the  middle  panel  is  the  calculated  stapes  displacement,  and 
the  lower  panel  is  the  envelope  of  basilar  membrane  displacements  and  the 
hazard  function  (explained  in  the  text).  The  arrow  in  the  upper  panel  indi¬ 
cates  the  peak  in  the  pressure  history  responsible  for  the  peak  BM  displace¬ 
ment. 
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FIG  1 1.  160-dB  howitzer  impulse  (Price.  1986a)  The  upper  panel  is  the 
pressure  history,  the  middle  panel  is  the  calculated  stapes  displacement,  and 
the  lower  panel  is  the  envelope  of  basilar  membrane  displacements  and  the 
hazard  function  (explained  in  the  text). 

feet  was  observed.  The  initial  peak  pressure  was  not  respon¬ 
sible  for  the  maximum  displacement  within  the  envelope 
(the  BM  only  went  a  little  over  half  the  envelope  height  in 
midcochlea),  rather  the  largest  displacement  was  a  result  of 
the  relatively  insignificant  peak  marked  with  an  arrow  in  the 
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FIG.  12.  154-dB  howitzer  impulse  (Price,  1986a).  The  upper  panel  is  the 
pressure  history,  the  middle  panel  is  the  calculated  stapes  displacement,  and 
the  lower  panel  is  the  envelope  of  basilar  membrane  displacements  and  the 
hazard  function  (explained  in  the  text). 
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FIG.  13.  1  JJ-dB  rifle  impulse  (Price,  1983b).  The  upper  panel  is  the  pres¬ 
sure  history,  the  middle  panel  is  the  calculated  stapes  displacement,  and  the 
lower  panel  is  the  envelope  of  basilar  membrane  displacements  and  the  haz- 
ard  function  (explained  in  the  text). 


FIG.  15.  Patterson  low-peak  impulse  (Patterson  el  at.,  1986).  The  upper 
panel  is  the  pressure  history,  the  middle  panel  is  the  calculated  stapes  dis¬ 
placement.  and  the  lower  panel  is  the  envelope  of  basilar  membrane  dis¬ 
placements  and  the  hazard  function  (explained  in  the  text). 


pressure  waveform.  We  presume  that  it  was  a  particularly 
effective  driver  because  it  occurred  near  ambient  pressure 
( when  the  stapes  was  not  clipping )  and  its  duration  was  such 
that  it  put  its  energy  where  the  ear  was  most  easily  driven. 


FIG.  14.  Patterson  hi-peak  impulse  (Patterson  etai,  1986).  The  upper  pan¬ 
el  is  the  pressure  history,  the  middle  panel  is  the  calculated  stapes  displace¬ 
ment,  and  the  lower  panel  is  the  envelope  of  basilar  membrane  displace¬ 
ments  and  the  hazard  function  (explained  in  the  text). 


The  evolution  of  the  hazard  function  was  also  instruc¬ 
tive.  First,  the  maximum  hazard  was  calculated  to  occur  in 
the  middle  of  the  cochlea.  This  was  true  for  all  the  impulses 
tested,  regardless  of  their  spectral  peak,  and  was  apparently 
due  to  the  fact  that  the  largest  oscillations  often  occured 
there  as  well  as  a  greater  number  of  the  smaller  oscillations. 
The  apical  portion  of  the  BM  would  typically  go  through  a 
large  single  displacement  when  significant  low-frequency 
energy  was  present  (reflected  in  the  envelope  of  displace¬ 
ment).  However,  this  apical  displacement  added  less  than 
might  be  expected  to  the  hazard  calculation  because  it  oc¬ 
curred  less  often  and  the  hazard  calculation  weighted  it  less 
because  it  occurred  in  a  wide  portion  of  the  BM  and  there¬ 
fore  it  represented  a  smaller  strain.  The  peak  of  the  hazard 
calculation  was  in  the  middle  of  the  cochlea  and  for  this 
impulse  reached  more  than  60,000  “units,”  mostly  during 
the  latter  portions  of  the  impulse  (after  10  or  15  ms). 

The  foregoing  calculation  compares  favorably  in  several 
respects  with  the  measured  hazard  (Price  et  al.,  1989).  Se¬ 
vere  losses  did  result  from  exposure  to  this  impulse,  and  the 
electrophysiological  and  histological  loss  measures  indicat¬ 
ed  damage  primarily  in  the  middle  of  the  cochlea. 

2.  The  160-dB  howitzer  Impulse 

Figure  1 1  is  a  parallel  to  the  previous  figure.  The  clip¬ 
ping  is  a  little  less  visible  in  the  stapes  displacement  as  would 
be  expected,  given  the  lower  peak  pressure  of  this  impulse. 
The  displacement  envelope  for  the  BM  shows  a  large  peak 
apically  and  a  smaller  one  in  midcochlea.  The  hazard  func¬ 
tion,  on  the  other  hand,  shows  the  reverse.  The  reason  for 
this  is  apparent  when  the  evolution  of  the  pattern  is  ob- 
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served.  The  apical  regions  of  the  cochlea  oscillate  less  often 
and  the  oscillations  are  weighted  less.  The  hazard  function 
grew  to  about  half  scale  based  on  the  response  to  the  first  5- 
10  ms  of  the  impulse  and  the  remainder  of  the  wave  roughly 
doubled  the  height. 

Again,  certain  of  the  model  calculations  match  the  loss 
data.  This  impulse  was  indeed  less  hazardous  than  the  166- 
dB  impulse,  and  losses  tended  to  occur  in  midcochlea. 

3.  The  154-dB  howitzer  Impulse 

The  model’s  response  to  the  154-dB  howitzer  impulse  is 
portrayed  in  Fig.  12,  where  we  see  a  continuation  of  the 
trends  established  in  the  two  higher  level  impulses.  For  this 
impulse,  the  hazard  function  reached  essentially  its  full 
height  after  the  first  phase  of  the  impulse  had  passed,  the 
latter  portions  of  the  wave  contributing  little  to  its  growth. 

The  match  between  the  hearing  loss  data  and  the  hazard 
calculation  for  this  impulse  raises  questions  of  a  more  gen¬ 
eral  sort.  The  greatest  losses  were  in  fact  in  the  midrange, 
which  does  match  the  hazard  calculation.  However,  this  im¬ 
pulse  and  the  160-dB  impulse  produced  essentially  the  same 
losses  in  the  cat  ear.  That  in  itself  is  somewhat  surprising 
[noted  in  the  original  article  ( Price  et  at.,  1989)  | ;  however, 
if  the  hazard  calculation  is  examined  after  the  first  10  ms, 
then  the  two  impulses  are  rated  about  equal.  We  tentatively 
speculate  that  if  the  middle  ear  muscles  begin  to  act  at  about 
this  point,  then  the  hazard  may  be  reduced  for  the  latter 
portions  of  the  wave.  This,  of  course,  would  be  true  for  the 
166-dB  impulse  as  well. 

4.  155-dB  rifle  Impulse 

In  Fig.  13,  the  data  for  a  rifle  impulse  are  depicted.  The 
pressure  and  time  scales  were  expanded  to  show  the  wave¬ 
form  better;  but  the  cochlear  scales  were  unchanged.  Some 
clipping  can  be  seen  in  the  stapes  displacement.  The  enve¬ 
lope  of  BM  displacement  shows  less  activity  at  the  apex,  as 
one  would  expect,  and  the  hazard  function  shows  its  peak  in 
midcochlea.  The  rifle  impulse  is  short,  therefore  the  hazard 
calculation  grows  only  in  the  first  few  milliseconds. 

The  match  between  calculated  and  measured  hazard 
(Price,  1983b)  is  relatively  good.  The  losses  were  in  fact 
primarily  in  the  middle  of  the  cochlea,  and  this  impulse  was 
about  as  hazardous  as  the  1 60-dB  cannon  impulse.  If  we  only 
consider  the  first  5-10  ms  of  the  cannon  impulse,  then  this 
relationship  also  holds  for  the  calculated  hazard.  If  the  en¬ 
tire  waveform  is  considered,  then  the  160-dB  howitzer  im¬ 
pulse  would  be  rated  as  much  more  hazardous  by  the  model. 

5.  Patterson's  hi- peak  and  low-peak  Impulses 

Patterson  and  his  colleagues  ( 1986)  did  an  interesting 
experiment  in  which  two  impulses  were  used  for  a  noise  ex¬ 
posure.  The  impulses  were  digitally  processed  so  that  at  the 
acoustic  output  of  a  speaker,  the  impulses  had  identical  mag¬ 
nitude  spectra  (peak  energy  at  1.4  kHz);  but  because  the 
phase  relationships  of  the  components  were  different,  the 
peak  pressures  differed  by  8  dB  (producing  wfiat  was  re¬ 
ferred  to  as  a  hi-peak  and  a  low-peak  impulse).  Both  histo¬ 
logical  analyses  and  behaviorally  measured  audiograms 


were  available  for  these  ears.  These  data  provided  an  inter¬ 
esting  challenge  for  the  model;  although,  strictly  speaking, 
the  data  were  gathered  on  the  chinchilla,  and  the  model  val¬ 
ues  were  for  the  cat. 

The  model’s  response  to  the  hi-peak  impulse  at  147  dB  is 
shown  in  Fig.  14.  At  this  pressure,  there  is  little  clipping.  The 
envelope  of  BM  displacement  and  the  hazard  function  both 
show  maxima  in  midcochlea,  as  one  would  expect.  The  mod¬ 
el’s  response  to  the  low-peak  impulse  at  139  dB  is  presented 
in  Fig.  15.  The  maximums  of  the  envelope  of  the  BM  dis¬ 
placement  and  the  hazard  function  are  both  somewhat 
smaller  than  for  the  high-peak  impulse  and  are  located  some¬ 
what  more  apically.  This  corresponds  with  their  effects  on 
the  ear.  The  hi-peak  impulse  produced  somewhat  more  PTS 
and  hair  cell  loss,  and  the  histological  damage  for  the  low- 
peak  impulse  was  located  somewhat  more  apically  than  the 
damage  for  the  high-peaked  impulse  (their  Figs.  15  and  16). 
Admittedly,  these  relatively  subtle  correspondences  may  be 
fortuitous;  but  we  note  them  with  interest. 

III.  INSIGHTS 
A.  Function  of  the  ear 

The  observation  that  the  ear  is  spectrally  tuned,  primar¬ 
ily  because  of  the  external  and  middle  ears,  is  hardly  new. 
Nevertheless,  it  is  a  concept  that  has  utility  in  explaining 
effect  of  noise  on  the  ear.  Certainly,  the  fact  that  the  middle 
of  the  cochlea  is  most  susceptible  to  damage  is  largely  due  to 
the  fact  that  the  transmission  path  works  most  efficiently  in 
the  midrange  of  frequencies. 

The  effect  of  a  nonlinear  stapes  has  not  been  so  thor¬ 
oughly  explored,  however.  A  limit  to  displacement  does  ex¬ 
plain  why  it  is  that  damage  is  less  than  we  might  otherwise 
expect  from  intense  impulses.  At  a  somewhat  more  subtle 
level,  a  displacement  limit  means  that  we  need  to  reexamine 
traditional  methods  of  rating  waveforms  for  hazard.  The  ini¬ 
tial  positive  peak  has  received  almost  all  the  attention  in 
damage  risk  criteria  (DRCs)  (it  contains  most  of  the  ener¬ 
gy  ) ;  yet  once  the  stapes  has  reached  its  limit,  the  ear  is  func¬ 
tionally  disconnected  from  the  acoustic  field.  Thus  the  fea¬ 
ture  that  controls  the  rating  of  hazard  may  be  a  minor  part  of 
the  stimulus.  On  the  other  hand,  when  the  acoustic  pressure 
is  near  ambient,  the  stapes  can  drive  the  ear  very  effectively. 
This  means  that  the  apparently  insignificant  “tails”  of  in¬ 
tense  impulses  may  actually  be  the  real  sources  of  cochlear 
input.  Furthermore,  when  the  initial  positive  pressure  ar¬ 
rives,  the  stapes  usually  starts  from  a  rest  position  and  can 
displace  20  //m  or  so.  However,  it  is  capable  of  displacing 
twice  that  distance,  if  driven  from  extreme  to  extreme,  which 
presumably  would  be  a  worst  case  condition. 

Although  not  presently  a  part  of  the  model,  the  middle 
ear  muscles  also  act  to  change  the  transmission  path.  They 
are  often  v  iewed  as  contracting  refiexively  to  loud  sound  and 
hence,  for  impulses  lasting  only  a  few  milliseconds,  as  being  a 
negligible  influence.  Alternatively,  given  their  short  latency 
of  contraction,  they  may  be  able  to  influence  the  later  por¬ 
tions  of  longer  impulses,  even  if  they  only  contract  reflexive- 
ly.  However,  they  may  be  active  for  other  reasons  and  their 
influence  should  perhaps  not  be  ruled  out  for  any  noise  expo¬ 
sure. 
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B.  Ameliorating  hazard 

Several  of  the  interesting  implications  from  these  calcu¬ 
lations  relate  to  novel  methods  of  ameliorating  hazard.  The 
presence  of  low  frequencies  at  the  input  to  the  ear  might  be 
more  of  a  help  than  a  hazard,  if  they  push  the  stapes  into 
clipping.  In  fact,  enough  pressure  in  the  outer  ear  or  middle 
ear  could  do  the  same  thing.  Small  barriers  between  the  im¬ 
pulse  source  and  the  ear  could  create  a  sound  shadow,  espe¬ 
cially  for  the  mid  and  high  frequencies  and  offer  significant 
protection  to  the  ear.  In  any  event,  because  the  model  is 
theoretically  based,  it  does  offer  the  possibility  for  evaluating 
and  developing  schemes  for  reducing  hazard. 

C.  Methods  of  assessing  hazard 

The  foregoing  analysis  suggests  that  because  of  the  non¬ 
linearity  in  the  middle  ear,  virtually  all  conventional  ap¬ 
proaches  to  the  rating  of  hazard  from  intense  sounds  (above 
140  dB)  are  likely  to  have  difficulties  that  may  be  insuper¬ 
able.  On  the  other  hand,  we  are  encouraged  with  the  suc¬ 
cesses  the  model  has  had  and  the  insights  it  is  providing. 
Should  additional  development  of  the  model  and  the  hazard 
function  continue  to  match  hearing  losses  already  measured 
and  predict  hazard  for  different  situations,  then  it  could 
serve  as  a  basis  for  a  new  DRC  for  impulse  noise.  The  use  of  a 
calculational  procedure  such  as  this  in  a  standard  is  not  tra¬ 
ditional;  but  given  the  modest  computational  requirements 
and  the  pervasiveness  of  the  PC-level  machines,  it  may  suc¬ 
ceed.  Alternatively,  it  is  conceivable  that  a  circuit  or  chip 
could  be  integrated  into  a  “hazard  meter,”  much  like  the 
current  sound-level  meters. 

The  foregoing  discussion  has  shown  briefly  that  the 
mathematical  models  are  reproducing  the  essential  elements 
of  hearing  hazard  and  arc  doing  it  in  a  theoretically  consis¬ 
tent  manner.  At  this  time,  no  other  method  of  rating  hazard 
comes  close  to  being  accurate.  We  are  very  much  encour¬ 
aged  that  these  models  are  on  the  right  track  and  will  prove 
to  be  sound  basis  for  a  new  design  criterion  for  impulse  pro¬ 
ducing  weapons. 
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